function [ylimvec]=limsaj(xmat,offset,hardbound) 
ymin=min(min(xmat));
ymax=max(max(xmat));
ylimvec=zeros(2,1); 
if nargin < 2 || isempty(offset)
    offset=0; 
end 
ylimvec(1)=ymin*(1-sign(ymin)*offset); 
ylimvec(2)=ymax*(1+sign(ymax)*offset); 
if nargin > 2 &&  ~isempty( hardbound ); 
    ylimvec(1)=min( [min(ylimvec(1));hardbound(1)] );
    ylimvec(2)=max( [max(ylimvec(2));hardbound(2)] );
end 
